Majorana single-charge transistor 



R. Hiitzen, 1 A. Zazunov, 1 B. Braunecker, 2 A. Levy Yeyati, 2 and R. Egger 1 

1 Institut fur Theoretische Physik, Heinrich- Heine- Universitat, D-40225 Dusseldorf, Germany 
2 Departamento de Fisica Teorica de la Materia Condensada C-V and Instituto Nicolas Cabrera, 
Universidad Autonoma de Madrid, E-28049 Madrid, Spam 
(Dated: March 8, 2013) 

We study transport through a Coulomb blockaded topologically nontrivial superconducting wire 
(with Majorana end states) contacted by metallic leads. An exact formula for the current through 
this interacting Majorana single-charge transistor is derived in terms of wire spectral functions. A 
comprehensive picture follows from three different approaches. We find Coulomb oscillations with 
universal halving of the finite-temperature peak conductance under strong blockade conditions, 
where the valley conductance mainly comes from elastic cotunneling. The nonlinear conductance 
exhibits finite-voltage sidebands due to anomalous tunneling involving Cooper pair splitting. 



PACS numbers: 71.10.Pm, 73.23.-b, 74.50.+r 

Introduction. — Topologically nontrivial insulators and 
superconductors exhibit many remarkable non-local fea- 
tures such as teleportation or non-Abelian statistics [1 2J. 
For a one-dimensional topological superconductor (TS) 
wire, such effects can be traced back to the existence 
of a zero-energy Majorana bound state (MBS) localized 
at each end [3HZ]- When a grounded TS is weakly con- 
tacted by a normal metal, the MBS is expected to pro- 
duce a characteristic zero-bias anomaly peak in the tun- 
nel conductance |8HT2]. Very recently, such a feature 
has been experimentally observed in tunnel spectroscopy 
using InSb or In As nanowires |T3l - fl"6] . where Majorana 
fermions are theoretically expected due to the interplay of 
strong spin-orbit coupling, Zeeman field, and proximity- 
induced superconducting pairing |17H19j . Recent reviews 
[U 121 EE US] have also summarized alternative MBS pro- 
posals. Here we discuss an interacting variant of previ- 
ously studied Majorana wire set-ups, the floating "Ma- 
jorana single-charge transistor" (MSCT) schematically 
shown in Fig. [T] A comprehensive picture of its trans- 
port properties in the presence of interactions emerges 
from our analysis below. Noting that the experimentally 
observed peak features could be related to a disorder- 
induced spectral peak |201I21| . our results should help to 
distinguish the Majorana state from alternative explana- 
tions in future experiments. 

Previous works [22-24 have studied electron-electron 
interactions in an isolated TS wire and found that Ma- 
joranas still exist under rather general conditions. As 
sketched in Fig. [T] we instead study a generalization of 
the set-up in Ref. |13| . where source and drain metal- 
lic electrodes contact the TS wire. We stress that the 
MSCT could be realized not only with nanowires but 
using most other Majorana proposals as well. In such 
a geometry, Coulomb blockade effects due to the finite 
charging energy E c of the TS can play a decisive role. 
For instance, one expects Coulomb oscillations of the con- 
ductance as a function of a gate voltage parameter n g , 
with peaks (valleys) near half-integer (integer) n g , while 



in the noninteracting (E c — 0) limit, the MBSs pinned to 
zero energy cause resonant Andreev reflection (AR) [8- 
lllj . with n g -independent linear conductance G = 2e 2 /h 
at temperature T — 0. Resonant AR also survives for 
E c < r = Tt, + albeit with reduced conductance |25j . 
For E c ^> r, Coulomb blockade is firmly established, and 
the peak conductance approaches the (spinless) resonant 
tunneling value G = e 2 /h, which has been pointed out 
as a signature of electron teleportation [SB] . 

In this paper, we consider Coulomb blockaded charge 
transport through the MSCT; for a variant with one su- 
perconducting and one metallic lead, see Ref. [27|. We 
provide an exact expression for the current in this inter- 
acting system, and develop three different approximation 
schemes to study Coulomb oscillations in the MSCT both 
for T — and finite T. We quantitatively describe the 
T = crossover of the peak conductance from G = 2e 2 /h 
to e 2 /h as E c /T increases, which constitutes a charac- 
teristic signature of Majoranas. Remarkably, this "halv- 
ing" of the peak conductance is universal and found to 
hold for arbitrary T. For the valley conductance, we find 
that elastic cotunneling dominates while AR is sublead- 
ing. We predict finite-voltage sidebands in the nonlin- 
ear differential conductance which are directly related to 
anomalous tunneling processes where the Majorana state 
and the Cooper pair number change simultaneously. The 
presence of Majoranas can be unambiguously identified 
in experiments by the magnetic field dependence of the 
sideband location. 

Model. —The MSCT Hamiltonian, H = H C + H t +H u 
contains a piece H c describing the TS wire, a tunnel 
Hamiltonian H t connecting the TS to the left (j = L) 
and right (j = R) electrode, and a term Hi describing 
the leads (we often use units with e = h = kg = 1). 
Topological arguments warrant that the TS wire holds 
a single unpaired MBS near each end [51 [7] described 
by the Majorana operator = 7J with anticommutator 
{ijiJj'} ~ We introduce the non-local fermion op- 
erator d = {jL+ijii)/ V% such that 7^ = (d+cf^/v^ and 
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Figure 1: (Color online) Majorana single-charge transistor 
(MSCT): The TS wire with Majorana end states is tunnel- 
coupled (Yl,Tr) to normal metal electrodes and Josephson 
coupled (Ej) to another bulk superconductor. Capacitive 
charging effects are encoded by E c and can be tuned by a 
gate voltage parameter n g oc V g . 



7j? = — i(d — S)/y2. With rid = cftd and the number op- 
erator N for Cooper pairs in the TS, the instantaneous 
charge state of the wire is described by (JV, rid), where 
the integer N and rid = 0,1 are eigenvalues of N and fid, 
respectively. With the phase \ conjugate to N, where 
[x, N] = i and e~ %x (e* x ) lowers (raises) N by one unit, 
we have 

H c = E C (2N + f\ d - n g f - Ej cos( X - </> s ). (1) 

The TS wire is assumed sufficiently long to exclude a 
direct tunnel coupling between 7^ and 7^. However, 
note that E c introduces a dynamical coupling between 
the Majoranas. Proximity-induced pairing correlations 
are required for MBS formation, and in Eq. ([I]) we in- 
clude Cooper pair exchange (with Josephson coupling 
Ej) between the TS condensate and another bulk su- 
perconductor (with fixed phase 4>s) |271 12"5] , We focus 
on the most interesting case of a large proximity gap 
Ats > max(i? c , T, T), where charge transport involves 
MBSs and the contribution of quasi-particles above the 
gap can be neglected. Next, electrons in lead j corre- 
spond to free fermions with chemical potential tij and 
(effectively spinless [25 J fermionic operators Cj } k for mo- 
mentum k. Hi is treated within the usual wide-band ap- 
proximation |29| and the bias voltage is eV = \ij, — [Ir. 
Taking into account charge conservation and expressing 
the Majoranas in terms of the non-local d fermion, the 
tunnel Hamiltonian reads 125 



r h = -L(d + S3 e-^), (2) 



where Cj — J2k c j,k, s l = +1 and sr — —1, and Xl,r de- 
notes the respective tunnel matrix elements Tunnel- 
ing from the TS to lead j thus proceeds either by destroy- 
ing the d state without changing N ("normal" tunneling) 
or by occupying the d state and simultaneously splitting 
a Cooper pair ("anomalous" tunneling), plus the conju- 
gate processes. Below we use the hybridization scales 



Tj = 2m>j\\j\ 2 , where Vj is the density of states in lead 
j. Experimentally, the Tj (and n g ) can be changed via 
gate voltages [T3] . 

Exact expression for current. — Using non-equilibrium 
Green's function (GF) techniques |50l I3T] . the current 
Ij flowing from lead j to the TS can be expressed in 
terms of the Keldysh GF G v .(t,t') = -i{Tcr]j (*)%(*'))> 
where 7c denotes Keldysh time ordering and the pseudo- 
fermion rjj has been defined in Eq. ^ . With the Fourier- 
transformed retarded, G^. (e), and Keldysh, G^. (e), com- 
ponents of G v ., we obtain Ij — (eTj/h) J de[F(e — 
»j) ImG*(e) + (i/2)G*(e)], where F(e) = 1 - 2/ = 
tanh(e/2T) encodes the Fermi function /(e) in the leads. 
Next we note that G^ (t, t) = as a consequence of 
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1/2. Hence we find the exact result 



j deF{e- H ) ImG«(e), 



(3) 



stating that the current can be computed from the spec- 
tral function oc ImG^ . The well-known expression for 
interacting quantum dots [32J has thereby been extended 
to the interacting Majorana wire; note that there are 
two spectral functions associated to the currents Ij, and 
In. Current conservation here implies Ij, + Ir + Is = 0, 
with the supercurrent 1$ flowing through the interface 
to the bulk superconductor. Below we define the con- 
ductance G = dl/dV using the symmetrized current 
I = [Ijj — Ir)/2. For E c = 0, the pseudo-fermions rjj re- 
duce to Majorana fermions jj , and the Lorentzian spec- 
tral function, — ImG^.(e) = Tj/(e 2 + T 2 ), implies reso- 
nant AR with G = 2e 2 /h [SHIP]. For finite E c , we shall 
present several complementary approximations in order 
to achieve a broad physical understanding of the MSCT 
transport properties. Equation ^ should also be useful 
for numerically exact calculations, e.g., using the numer- 
ical or density-matrix renormalization group. 

Equation- of-motion (EOM) approach. — We con- 
structed an EOM approach for G^. to access 
the linear conductance near a peak. Within 
this method, we introduce the Nambu spinors 
^d = (d, e~' x tff) and the corresponding retarded 



-*e(* - t)({* d (t),%{t)}). The EOM 



GF, G« : 

for G dd generates higher-order GFs of the type 
t')({N m (t)«£ d (t),*l(t')}), which 



1 N™dd 



-iQ(t 

we truncate at the level m — 2 and solve in a self- 
consistent way [33]. The resulting GF G dd then yields 
G^. = |Tr [(1 + Sj(T x ) Gf d ] with Pauli matrix a x . 
Finally, we obtain the conductance from Eq. Q3J) . This 
approximation is valid by construction for E c > T, but 
the imposed self-consistency [33] allows us to extend 
it to E c < r, where the resulting conductance (being 
determined by trunctated fluctuations) gives a lower 
bound for the exact result. 

Zero-bandwidth model (ZBWM). — Next we study the 
ZBWM where each lead is represented by just a single 
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fcrmion site and only a finite number of TS Cooper pairs 
(N < iV max ) is included. The Hilbert space then has 
the finite dimension 8iV max , which allows us to numer- 
ically calculate the spectral density oc ImG,^ (e) via its 
Lehmann representation, with poles phenomenologically 
broadened by T. A similar description has been pursued 
before for E c = Ej = [34 . With this spectral function, 
Eq. Q yields the conductance within the ZBWM. 

Master equation and cotunneling processes. — For T > 
r, the GF formulation reduces to a master equation de- 
scription including sequential tunneling and cotunneling 
processes (for simplicity, Ej = here). The station- 
ary probability distribution Pq for having Q = 2N + n<j 
particles on the TS then obeys J2q'^q[Pq'Wq'^q — 
PqWq-^q'] = 0. All non-vanishing transition rates 
Wq->.qi are specified in terms of rates obtained under 
a systematic second-order T-matrix expansion in Tl^r 
|35| . With the electrostatic energy Eq = E C (Q — 

n g ) 2 , sequential tunneling yields the rate ^jQl^Q-i-i = 
(Tj/2)f(EQ±i — Eq =p fij) for one particle tunneling into 
(out of) the TS from (to) lead j = L, R. Next, elastic co- 
tunneling transfers a particle from lead j to the opposite 
lead —j with virtual excitation of the TS states Q ± 1. 
The elastic cotunneling rate is 

lf Q C) = ^ / def(e - - f(e - (4) 



e-(E Q+1 -E Q )+iO e - (E Q - Eq^) - iO 

where the two terms come from the interference of normal 
and anomalous tunneling. We note in passing that for 
large Ats> inelastic cotunneling does not contribute at 
all, while the conventional elastic cotunneling rate due to 
quasi-particle states above the gap (and without MBSs) 
would be much smaller, r( EC ) cx T l T r /A ts [3D]. To the 
same order in Tl,,r, we also have local (and crossed) AR 
processes, where an electron and a hole from the same 
(different) lead(s) are combined to form a Cooper pair, 
Q —> Q + 2; the reverse process describes Cooper pair 
splitting, Some algebra yields the AR rates 



_ 1 + h-3' T 3 T 3' 



r (AR) 



de / de' 



(5) 



x/(±(e - H ))f{±{e' - ^))S(e + e' T (E Q±2 - E Q )) 



s j s j' 



ct(Eq ±1 -Eq) + iO e> t (Eq±i - Eq) + iO 

where j — j' (j ^ j') corresponds to local (crossed) AR. 
The iO terms indicate that regularization of the integrals 
in Eqs. Q and (|5| is necessary. Application of the gen- 
eral regularization scheme in Refs. (313 EZ] then implies 
that the principal value of these integrals needs to be 
taken. Given these rates and the (numerical) solution 
Pq of the master equation, the currents Ij then readily 
follow. 
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Figure 2: (Color online) Coulomb oscillations in the MSCT. 
Main panel: Conductance G vs n g from master equation for 
Ej = 0, T = 2r and several E c . Inset: Same but from ZBWM 
for T = 0, E c = 5r and several Ej. 




Figure 3: (Color online) Peak conductance G vs E c /T on 
a semi-logarithmic scale. Main panel: Comparison of T = 
results using perturbation theory in E c /F [25\ (blue solid 
curve), the EOM approach (red dotted-solid curve), and the 
ZBWM (black dashed curve). The shown EOM results are 
quantitatively valid only for E c > V (solid part) but give a 
lower bound when E c < F (dotted part). Here Ej — since 
G only weakly depends on Ej for Ej < E c . Inset: Same but 
from master equation for several temperatures T > F 



Coulomb oscillations. — Let us first address the rig- 
dependence of the linear (V — > 0) conductance, see Fig. [2] 
we take = Tr = Y/2 in all figures. Both the master 
equation (main panel, finite T) and the ZBWM (inset, 
T = 0) reveal clear conductance oscillations in the MSCT 
for E c ^> r, with peaks (valleys) for half-integer (integer) 
gate voltage parameter n g . The main panel shows that 
the peak (valley) conductance is halved (strongly sup- 
pressed) when going from the noninteracting to the deep 
Coulomb blockade limit. For Ej = and (T,T) < E c , 
the lineshape of the valley conductance is obtained in 
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closed form, 
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(6) 



where 6 = n g — [n g ] with \6\ <C 1 is the deviation from 
a valley center. Equation ([6| comes from elastic cotun- 
neling, with constructive interference of the normal and 
anomalous tunneling contributions [see Eq. Q], while 
AR is strongly suppressed in this limit. The inset of Fig. [2] 
shows that G increases when the Josephson coupling 
Ej grows. One can understand this by noting that for 
Ej ^> E c , one ultimately reaches the resonant AR limit 
of a grounded TS, with the n g -independent T — con- 
ductance G = 2e 2 /h. We find that AR yields significant 
conductance contributions for Ej > E c . which are best 
detected through the non-local conductance dlhjd^R. 
However, we will discuss this quantity in detail elsewhere. 

Peak conductance. — Results for the peak conductance 
are shown in Fig. [3] For T = (main panel), we obtain 
the full crossover from G — 2e 2 /h to G — e 2 /h as E c /T 
increases. The known small--E c behavior [55] is nicely 
reproduced by the ZBWM calculation. In the opposite 
large- E c limit, the EOM method is very accurate and 
Fig. [3] suggests that the simple ZBWM already captures 
the crossover from resonant AR [SHIP] to electron tele- 
portation |26] surprisingly well. The inset of Fig. [3] again 
demonstrates the universal halving of the finite-T peak 
conductance, see also Fig. [2] Since experiments so far 
were conducted in the high-temperature regime T > T 
|13| . let us now specify the lineshape near a conductance 
peak for E c > T. Using S = n g - [n g ] - 1/2 with \S\ < 1 
for the deviation from a peak center, truncation of the 
master equation to two charge states gives 
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(7) 



We stress that the peak conductance [G pea k(<5 = 0)] is 
indeed halved compared to E c = |8J. Moreover, it 
exhibits both a thermal and an interaction-induced re- 
duction. 

Finite-voltage sidebands. — Next we discuss the differ- 
ential conductance at finite bias voltage V. Master equa- 
tion results for T = 2T are shown in Fig. |4j Start- 
ing with the main panel, we find sideband peaks when 
eV is equal to an integer multiple of AE C . For these 
voltages, the chemical potentials hl,r are resonant with 
two (almost) degenerate higher-order charge states, im- 
plying additional sequential tunneling contributions be- 
yond the resonant transition determining the linear con- 
ductance peak [Eq. Q]. Note that the fluctuations in 
N needed to reach higher-order charge states can only 
be achieved through anomalous tunneling processes [see 
Eq. ([2])]. Similar sideband peaks are also found for other 
n„; the integer- 7i„ valley case is shown in the inset of 
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Figure 4: (Color online) G = dl /dV vs voltage V from the 
master equation for T = 2r, Ej = 0, and several E c /Y. The 
main panel (inset) is for half-integer (integer) n g . 
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Figure 5: (Color online) Same as Fig. [4] but from ZBWM for 
T = 0, E c = 5r and several Ej. 



peaks as Ej is changed, determined from the ZBWM at 
T = 0. For half-integer n g , the sideband peak position 
observed in the main panel of Fig. [5] is well described by 
eV ~ 4^ c ^l + (Ej/2E C ) 2 , which comes from Josephson 
coupling between the two relevant charge states. Since 
Ej can be tuned by applying a small magnetic field par- 
allel to the junction between the TS and the bulk super- 
conductor, an experimental observation of the sideband 
peak and its shift with magnetic field [cf. the expression 
for the peak position above] would provide clear evidence 
for the anomalous tunneling mechanism, and thereby for 
the presence of MBSs. 

Conclusions. — In this paper, we have studied the 
transport properties of an interacting Majorana single- 
charge transistor. Our results should be directly rele- 
vant for experiments extending existing work, see, e.g., 
Ref. [13] , where conductance peaks for tunneling into Ma- 
jorana wires were reported. When a gate voltage param- 
eter n g is varied, we find Coulomb oscillations, where the 



Fig. |4] In Fig. [5] we show the evolution of the sideband behavior of the peak and valley conductance has been 
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characterized in detail. The Majorana fermions in this 
system could be identified by observing sideband peaks 
in the nonlinear conductance and from the crossover be- 
havior of the Coulomb peak conductance. — This work 
was supported by the DFG (Grant No. EG-96/9-1 and 
SFB TR 12), by the EU network SE2ND, and by the 
Spanish MICINN under contract FIS2008-04209. Note 
added: After submission of this work, we learned of an 
unpublished study of the MSCT model by L. Fu and C.L. 
Kane. 
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Supplementary material: EOM approach 



Using a notation as 



the main text, the 



current follows from the retarded GF G dd (t) — 
-i6(t)({* d (t), which is a 2x2 matrix in Nambu 

space with spinors ^>d = (d, e~ tx <v) T . In energy (e) 
space, it obeys the EOM 



(e - E a + V g + it )G& = 1 + AE C T R 

f E r 



(8) 



where V g = 2E c n g , E = ^ J , and T = J2j ^(1 + 

Sjcr x ). r^ dd in Eq. (|8| is the first in a hierarchy of vertex 
functions (to = 1,2,...), 



r 



N m dd 



( t ) = -iemN m m d (t),nm), (9) 



which are generated sequentially through their EOM. In 
particular, 

(e - E + V g + if)r* dd = A- iTG dd + 4E c T* 2dd , (10) 



with r = Y] j Tj Sj <r x and 







A= ({JV* d ,*U) = ( vv " I . (11) 

u dl dil \ (N - (1 - n d ))J 1 ; 

Next we discuss an approximation closing the above set of 
equations. We see from Eqs. ([8| and ( 10 1 that the energy 
dependence of each higher-order vertex function produces 
a pole in G dd , and the scale of the energy spacing between 
the poles is set by the prefactor AE C with which the vertex 
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functions appear in the EOMs. For E c > T, the physics 
is therefore controlled by a small number of poles only, 
which allows us to close the EOM hierarchy by truncation 
(effectively keeping just a few poles) . An approximation 
keeping only two poles can be achieved by imposing the 
variational condition 



r 1 * — HrK 

1 N 2 dd — nL Ndd' 



(12) 



with a Nambu matrix B. From this we obtain the closed 
equation 



(e - E + if - 4E C B + V g )(e - E + if + V g ) 



4iE r T 



G* i = e-E Q +it + 4E c (A-B) + V g . (13) 



For symmetric contacts, = T^, we have r — 0, f = T, 
all matrices become diagonal, and 



B- X A 



l-B~ l A 



dd 



e-E +ir- AE C B + V g e-E + iT + V g 

(14) 



The vertex function correspondingly reads 



1 Ndd 



A 



e - E Q + iT - AE C B + V 



(15) 



This EOM implementation is valid for gate voltages V g 
close to resonance for either the [G dd ]n or [G dd ]22 ma- 
trix entries. The on-resonance entry will then be dom- 
inated by a single central pole, while the other entry is 
off-resonant due to the 2E C energy shift between the two 
entries in the matrix Eq- Through this shift, the Fermi 
surface lies almost in the center between the two poles 
of G dd , and keeping just these two poles is sufficient for 
E c > r. For E c < r, however, higher-order poles become 



important. These come from higher-order fluctuations of 
N m in the neglected vertex functions. Hence Eq. ( 14 I can 



be interpreted as a truncation of fluctuations in the num- 
ber of Cooper pairs. Therefore, while we obtain quanti- 
tatively accurate conductance results for E c > T, due to 
the truncation of Cooper pair fluctuations we get only a 
lower bound for the conductance when E c < T. Within 
the restrictions imposed by the truncation, however, we 
achieve an optimal solution for G dd by exact fulfillment 
of the following sum rules. 

From Eqs. (14) and (15 1, a self-consistent computa- 



tion of four parameters is required, namely (N) and (n d ) 
appearing in the matrix A [see Eq. (11)], and the two 
diagonal matrix entries B\\ and -B22 of B [see Eq. ( 12 1] . 



(Alternatively, it can be advantageous to fix (rid) but ad- 
just V g self-consistently. ) These values are determined 
from exact sum rules, 



def d (e)lm[Gf d (e)}u = (h d ), 



- / def d (e)lm[Gf d (e)] 22 = 1 - (n d ), 



(16) 
(17) 



de/ d (e)Tr{Imr^ d (e)} = (N - (1 - n d )), (18) 



where f d is the distribution function on the TS wire; in 
equilibrium, f d = j ' . (Note that through the trace over 
Imrj^^, averages of the form (Nn d ) cancel out.) How- 



ever, when using Eq. (14 1, we have to self-consistently 
adjust four parameters with three sum rules only. We 
thus impose TrB = B n + B 22 = (N) + (1 - (h d )), which 
reproduces in the large-_E c limit the resonances found 
with the ZBWM and master equation approaches around 
(n d ) = 1/2. The EOM results shown in Fig.[3]thcn follow 
(with T = and given ratio E c /T) by self-consistently 



solving for G 



